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Abstract: We herein utilize the general three-body problem (GTBP) as a model, in order to simulate 
resonant systems consisting of a star and two planets, where at least one of them is highly eccentric. 
We study them in terms of their long-term stability, via the construction of maps of dynamical stability 
and the computation of the corresponding families of periodic orbits. We identify the way their survival 
is connected with the regions of regular motion in phase space, which, in turn, were created by stable 
resonant periodic orbits in their vicinity. Gonsequently, a phase protection mechanism is provided and 
the planets avoid close encounters and collisions even on long timescales. We apply our methodology 
to the extrasolar system HD 82943. 

1 Introduction 

The number of known extrasolar planets has increased considerably over the last years and the dis¬ 
covery of these systems has raised many interesting questions e.g. about their formation, composition, 
habitability etc. Today, about 500 systems seem to consist of more than two planets and in many 
of them, the gravitational interaction between the planets is not negligible. Therefore, interesting 
questions about their dynamical evolution and stability are raised, too. 




(a) (b) 

Eigure 1: Distribution of eccentricities of exoplanets (data from exoplanet, eu) 

The long-term stability of multi-planet systems becomes very interesting when eccentricities are 
quite high. In this case, the almost Keplerian elliptical orbits can be close to each other or intersect 
and, therefore, mutual gravitational interactions may be very strong. A mechanism that prevents 
the close encounter of planets is the resonant phase protection. When planets revolve so that their 
mean motion ratio is almost rational ni/n 2 ~ p/q, P^q ^ N, i.e. we have a mean motion resonance, 
we can determine initial conditions of the planets, for which, although the orbits are very close to 
each other, the planets avoid coming close to one another and, therefore, their mutual interactions are 
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Figure 2: A dynamical map in the plane a — 62 for Jupiter-Saturn masses (ei = 0.3, zui = W 2 = 0, 
Ml = M 2 = 0, where Mi, M 2 are the mean anomalies). The corresponding resonances for the highly 
eccentric islands are indicated. Gray curve shows the collision line - above which the planetary orbits 
intersect. 


weak. Subsequently, planetary systems of high eccentric orbits can only survive if they are trapped in 
a resonance and have the appropriate phases. 

According to the data provided by the extrasolar planets encyclopedia, many exoplanets have large 
eccentricities. The distribution of the eccentricities with respect to the orbital period is given in Fig. 
I^a). An important question is whether such eccentric planets can have companions [7]. Nevertheless, 
according to the analysis of radial velocity data, multiple systems with eccentric planets seem to exist 
e.g. the systems HD 82943, HD 3651, HD 7449, HD 107222 and HD 89744 (see Fig. [^b)). Data 
from Kepler also indicate many planets with high eccentric orbits [3]. In this paper, we present a 
methodology for obtaining stable orbits in two-planet systems where, at least, one of the planets 
revolves in a highly eccentric orbit. The main aspects of two-planet systems dynamics are presented 
in 13 . 

2 Methodology 

We model a two-planet system with the general planar three body problem consisting of a very massive 
star with mass rriQ and two planets with masses mi and m 2 (the inner and the outer planet, respec¬ 
tively). By using the osculating orbital elements, the parameters that define explicitly the dynamical 
position of the system are the ratio of semimajor axes a = a 2 /ai, the eccentricities ei and 62 , the angle 
of pericenters Atu = zj 2 — txji and the initial relative position on their osculating ellipses, which can 
be given by the resonant angle 0 = pX 2 — qXl — {p — q)'^i [I]- 

As we mentioned above, highly eccentric orbits of long-term stability can exist in resonant regions. 
The centers of such regions (exact resonances) are periodic orbits of the three body problem in a rotating 
frame of reference [ 2 ]. Resonant periodic orbits are not isolated, but they form monoparametric families 
in phase space. For the case of symmetric periodic orbits Atu and 0 are constant and take the values 0 
or TT. This means that the planetary orbits are aligned or anti-aligned and the planets are in conjaction. 
These families can be represented by characteristic curves in the plane ei — 62 . However, along such 
curves a varies slightly. 

With the term “stability” we refer to orbits that evolve regularly, namely they are quasiperiodic 
orbits winding invariant tori in phase space. Chaotic orbits may also survive for long time integrations, 
however, the majority of them terminates at a collision or escape of one of the planets. Chaos indicators 
can be used for detecting the chaotic nature of orbits in relatively short time intervals (e.g. see 0 ). 
In our study, we use the de-trented FLI defined as the function 

DFLI{t)=supt{m\/t), 

where ^ is the deviation vector of the orbit, which is evaluated along the numerical integration of 
the orbit simultaneously with the integration of the linearized system. Dynamical maps of stability 
can be computed by considering planar grids of initial conditions, numerical integration of orbits and 


12th Hel.A.S Conferenee 


2 












Antoniadou and Voyatzis: Highly eccentric exoplanets 


1 


0.5 


0 0 


- 0.5 


-1 

-1 - 0.5 0 0.5 1 

ei 

Figure 3: Dynamical map for the 2/1 resonant motion and for the masses of the planets HD 82943b,c. 
Each quadrant corresponds to the indicated configuration (^, Atu). Bold blue and red curves are the 
the stable and unstable families of periodic orbits, respectively. The gray bold line is the collision line. 

computation of the value DFLI = DFLI{tmax) El- An example is shown in Fig. where a map 
on the plane a — 62 is given by fixing ei = 0.3, Atu = 0 = 0° and setting the masses mi = 0.001, 
m 2 = 0.0003 (Jupiter - Saturn planets). Integrations performed for 200 Ky (with ai = lAU). Dark 
colors {DFLI < 10^) correspond to long-term stable orbits, while yellow color indicates strongly chaotic 
orbits with close encounters in relatively short time intervals. A notable main property of the DFLI is 
the sharp distinction between chaos and order. By obtaining a large value e.g. DFLI > 10^^ we can 
stop the integration and classify the orbit as chaotic. 

The dynamical map shows that for 62 < 0.4 (bellow the collision line) there exist many regions of 
stable motion. However, above the collision line, only distinct and small (but not negligible) stable 
regions are detected. These regions of initial conditions provide resonant motion and each one is built 
around a particular family of periodic orbits. 

3 An application to HD 82943 extrasolar system. 

The suggested parameters for the position of the system HD 82943 are [ 8 ] 

mi/mo = 0.0041, m 2 /mo = 0.0042, ai = 0.74, a 2 = 1.18, ei = 0.43, 62 = 0.2, Atu = 25°, 0 = —2°. 

The above values indicate that the system should be trapped in 2/1 resonance and close to the symmet¬ 
ric configuration Azu = 0 = 0°. For the particular masses we computed the main families of symmetric 
periodic orbits at the initial configurations ( 6 >, Atu) = (0,0), ( 0 , 7 r), ( 7 r, 0 ), ( 7 r, 7 r). For the same config¬ 
urations we computed the dynamical maps at the plane ei — 62 with a ~ 1.59. Each map provides a 
quadrant of the global map presented in Fig. The ei — 62 characteristic curves of the families of 
periodic orbits are also presented. Bold blue segments indicate stable periodic orbits, while red ones 
consist of unstable periodic orbits. It is clearly shown that the configurations (0, tt) and (tt, 0) do not 
support any stable orbits and families of stable periodic orbits do not exist. Instead, the presence of 
the stable family of periodic orbits for ( 6 >, Aw) = (tt, tt) is accompanied by a broad region of initial 
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conditions of regular orbits. Wholly this region is located above the collision line and stable orbits can 
be found even for ^ 1. 
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Figure 4: The evolution of orbital elements of the extrasolar system HD 82943. 

The position of HD 82943 in the map is indicated by the red-star like symbol. We obtain that the 
family of periodic orbits in the (0, 0) configuration is surrounded by a region of stability. Starting with 
initial conditions in this region, the evolution of orbital elements is regular, and, particularly, the angles 
0 and Azj librate around 0°. The system HD 82943 is located in this region of stability. Starting with 
the initial conditions, given above, the orbital elements evolve as it is shown in Fig. and this picture 
is expected to remain unaltered for many Gyears. 
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